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The  carbon  deposition  at  the  Triple  Phase  Boundary  (TPB)  of  the  Nickel/Yttrium-Stabilized  Zirconia  (YSZ) 
interface  is  studied  using  the  first-principles  method  based  on  density  functional  theory,  with  consid¬ 
eration  of  the  interface  oxygen  vacancy.  It  is  found  that  the  CH  fragment  (the  most  stable  dissociation 
products  of  CH4  on  Ni  catalyst)  can  easily  diffuse  and  be  trapped  at  the  O  vacancy.  The  trapped  CH  can 
dissociate  to  C  and  H  with  a  much  lower  dissociation  barrier  (0.74  eV)  as  compared  with  that  (1.39  eV)  on 
the  pure  Ni  (111)  surface.  Therefore,  we  propose  that  the  carbon  deposition  may  form  easily  at  the 
interface  oxygen  vacancy  of  TPB  as  compared  with  that  on  the  pure  Ni  (111)  surface,  which  offers  new 
understanding  on  the  carbon  deposition  of  the  Ni/YSZ  anode  of  solid  oxide  fuel  cell. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Solid  oxide  fuel  cells  (SOFCs)  are  expected  to  be  a  crucial  tech¬ 
nology  in  the  future  power  generation  [1,2].  SOFCs  offer  many 
desirable  advantages  compared  to  other  types  of  fuel  cells  and 
conversion  devices  due  to  the  use  of  solid  electrolytes,  lack  of 
moving  parts,  ability  to  circumvent  precious  metal  use,  high 
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efficiency,  low  pollution,  and  fuel  flexibility.  The  fuel  flexibility  is 
achieved  by  reducing  the  oxygen  molecules  in  the  cathode 
compartment  and  transporting  the  oxygen  ions  for  fuel  oxidation 
through  the  solid  state  electrolyte,  which  requires  high  operating 
temperatures  (600-1000  °C).  Unfortunately,  the  high  operating 
temperatures  of  SOFC  lead  to  a  number  of  complications  that  can  be 
detrimental  for  the  cell  performance.  These  include  low  anode 
stability  to  redox  cycling  [3,4],  adverse  reactions  between  adjacent 
cell  components  [5,6],  and  temperature  gradients  in  the  cell  parts 
occurring  during  the  start-up  and  shut-down  phase  [7,8].  More¬ 
over,  the  direct  electrochemical  oxidation  of  hydrocarbon  fuels  in 
SOFCs  is,  in  principle,  possible,  but  has  resulted  in  carbon 
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deposition  (coke  formation)  at  elevated  temperatures  (most 
severely  at  T  >  800  °C),  which  results  in  chemical  deactivation  of 
the  anode  9].  The  formation  of  coke  over  nickel-based  catalysts 
during  the  partial  oxidation  or  the  steam  and  CO2  reforming  of 
methane  or  other  hydrocarbons  has  been  extensively  investigated 
during  the  past  decades  [10-12].  Coke  accumulation  has  been 
suggested  to  be  a  balance  between  the  formation  and  the  removal 
of  coke  which  determines  the  degree  of  coking  [13  .  The  growth  of 
carbon  on  the  nickel  surface  during  a  hydrocarbon  reforming  pro¬ 
cess  has  been  explored  experimentally  in  detail  by  Sehested  [14]. 
Three  types  of  carbon  have  been  observed  in  a  reformer:  pyrolytic, 
encapsulating,  and  whisker  carbons.  Ni  catalyzes  the  formation  of 
carbon  deposition  at  the  anode  of  SOFCs  under  reducing  condition 
[15,16],  which  would  block  the  active  sites  and  deactivate  the  Ni 
catalysts,  ultimately  destroy  the  catalyst  completely  [16-20].  Coke 
formation  on  nickel  surfaces  has  been  fairly  well  investigated.  The 
carbon  growth  on  the  Ni  catalysts  is  also  generally  accepted  to 
occur  via  the  dissolution-precipitation  mechanism  21-23].  For 
example,  methane  (CH4)  adsorbs  onto  the  nickel  surface  and  then 
decomposes  into  C  and  H;  later,  the  carbon  atoms  dissolve  into  the 
bulk  of  the  Ni  particles,  diffuse  through  the  Ni  particles,  and  finally 
precipitate  as  carbon  at  their  rear  side.  Though  it  is  true  that  Ni 
generally  catalyzes  carbon  deposition  as  well  as  reforming  re¬ 
actions,  it  also  depends  on  many  different  parameters  such  as 
temperature,  Ni/YSZ  formulation,  use  of  additives  and  steam  to- 
carbon  ratio  24-26].  Steam  reforming  of  hydrocarbons  involves 
a  risk  of  carbon  deposition.  The  probability  for  carbon  depositions 
depends  on  the  steam-to-carbon  ratio  and  on  the  catalytic 
composition  [27,28].  The  risk  of  carbon  deposition  can  be  decreased 
if  the  pre-reforming  is  carried  out  before  the  cell,  or  the  anode 
structure  or  configuration  is  modified.  More  detailed  introductions 
to  the  carbon  deposition  studies  are  shown  in  several  excellent 
review  papers  [14,25,29-31]. 

In  general,  researchers  [32-34]  studied  the  carbon  formation 
started  at  the  dissociation  of  CFI4  on  the  Ni  (111)  surface.  The  suc¬ 
cessive  dehydrogenations  of  CFI4  to  CFI3,  CFI2,  CFI  are  favorable  on  Ni 
(111 )  with  the  activation  barriers  less  than  1  eV  (0.91, 0.67, 0.30  eV), 
while  the  barrier  of  last  step  of  dehydrogenation  of  CFI  is  as  large  as 
1.36  [32],  1.39  [33]  eV  with  endothermic  reaction  energies  of  0.54 
[34]  and  0.49  [32]  eV.  Most  of  the  theoretical  studies  of  methane 
adsorption  and  dissociation  on  Ni  catalysts  arrived  at  conclusions 
that  the  CH  species  has  the  most  thermodynamic  stability  [35-37] 
among  the  dissociation  species  (CHX  (x  =  0-3))  of  CFI4.  For  example, 
DFT  results  of  Nikolla  and  coworkers  [36]  on  the  reaction  energies 
for  various  elementary  steps  involved  in  methane  stream  reforming 
on  Ni(lll)  showed  that  the  CFI  fragment  is  the  most  adsorbate 
species  on  Ni(lll)  surface.  As  known,  the  anode  of  SOFC  is  the  Ni/ 
YSZ  composite,  instead  of  the  Ni  itself.  The  ideal  Ni(lll)  only  rep¬ 
resents  the  anode  region  beyond  the  TPB.  It  is  generally  believed 
that  the  active  sites  of  the  Ni/YSZ  anode  are  covered  with  deposited 
carbon,  resulting  in  the  deactivation,  loss  of  cell  performance  and 
lower  SOFC  reliability  [15,26,38,39].  Shishkin  and  Ziegler’s  studies 
[40,41  ]  on  hydrogen  oxidation  at  the  Ni/YSZ  interface  indicated  that 
the  most  active  sites  are  located  at  the  Ni/YSZ  interface.  Our  pre¬ 
vious  works  [42,43]  showed  that  the  sulfur  atom  would  be  strongly 
adsorbed  at  the  vacancy  site  at  the  Ni/YSZ  interface  and  is  more 
favorable  than  that  on  a  Ni  surface.  These  findings  all  indicate  that 
the  vacancy  site  at  the  Ni/YSZ  interface  is  highly  active  compared 
with  the  ideal  Ni  (111)  surface  beyond  the  TPB  region.  Recently, 
Shishkin  and  Ziegler  [44]  studied  the  carbon  removal  from  the 
anode  triple-phase  boundary  (TPB)  of  SOFCs  by  adsorbed  water 
molecules.  They  found  that  the  formation  of  interfacial  carbon  (at 
the  O  vacancy  site)  is  more  favorable  than  that  on  a  Ni  surface, 
indicating  that  the  TPB  of  Ni/YSZ  can  also  be  clogged  by  carbon 
atoms.  Ihara  et  al.  have  also  studied  [45]  the  deposited  pyrolytic 


carbon  at  the  three-phase  boundary  in  experiments.  However,  the 
studies  on  the  carbon  deposition  at  the  TPB  of  the  Ni/YSZ  interface, 
especially  for  that  with  O  vacancy,  are  very  scarce. 

In  this  work,  the  processes  of  CH  diffusion  to  the  vacancy  site  at 
the  Ni/YSZ  interface  and  the  dehydrogenation  of  CH  at  the  vacancy 
site  are  studied.  It  is  found  that  the  CH  can  dissociate  easily  at  the 
vacancy  site  with  a  much  lower  barrier  (0.74  eV)  as  compared  with 
that  (1.39  [33],  1.36  [32]  eV)  at  the  Ni  (111)  surface. 


2.  Model  and  computation  method 

All  calculations  presented  in  this  work  are  performed 
employing  the  periodic  density  functional  theory  (DFT)  method 
implemented  in  the  Vienna  Ab-Initio  Simulation  Package  (VASP) 
[46].  The  exchange-correlation  interactions  are  treated  with  the 
Perdew-Burke— Ernzerhof  (PBE)  functional  47].  Spin-polarized 
calculations  are  applied  throughout.  The  electron-ion  in¬ 
teractions  are  treated  using  the  projector  augmented  wave  (PAW) 
method  [48,49].  The  wave  functions  are  expanded  in  plane  waves 
with  a  cut  off  energy  of  408  eV.  The  TPB  model  of  the  Ni/YSZ 
cermet  used  in  the  Shishkin  and  Ziegler’s  work  [40]  is  adopted  as 
the  substrate,  where  the  YSZ  slab  has  the  horizontal  dimensions  of 
12.56  x  7.25  A.  In  the  proposed  structure,  both  Ni  and  YSZ  face 
each  other  by  the  (111)  crystallographic  planes,  with  a  small  lattice 
mismatch  of  3%  in  the  direction  with  sustained  translational 
symmetry.  Experimentally,  Abe  et  al.  [50]  and  other  researchers 
[51,52]  fabricated  and  characterized  the  Ni/YSZ  anode  cermet, 
which  has  the  structure  with  the  (111)  planes  of  the  Ni  part  par¬ 
allel  to  the  (111)  planes  of  YSZ.  Using  the  transmission  electron 
microscopy  technique  (TEM),  they  have  shown  a  clear  absence  of 
amorphous  phases  at  the  interface  with  a  (111 )//(lll )  orientation 
relationship  between  Ni  and  YSZ.  A  vacuum  layer  of  15  A  is  used  to 
separate  the  periodic  images  in  the  direction  perpendicular  to  the 
surface.  The  Monkhorst-Pack  k-point  mesh  of  2  x  3  x  1  is  used  for 
the  Brillouin  zone  (BZ)  sampling.  The  atoms  in  the  bottom 
multilayer  are  kept  fixed  for  all  calculations.  Structural  optimiza¬ 
tion  of  all  systems  is  performed  until  the  atomic  forces  drop  below 
0.02  eV  A-1.  The  climbing  image  nudged  elastic  band  (CI-NEB)  [53] 
method  is  employed  to  calculate  the  transition  states  and  migra¬ 
tion  barriers.  The  adsorption  energy  of  an  atom  or  molecular 
fragment  (e.g.,  C  or  CH)  is  defined  by 

^ads  =  ^CH  or  C  +  ^Ni/YSZ-Ov  “  ^CH  or  C-Ni/YSZ-Ov  (i) 

where  Ech  or  c  is  the  energy  of  a  single  CH  fragment  or  a  C  atom 
simulated  in  the  8  x  8  x  8  A  box;  ECh  or  c-Ni/Ysz-ov  and  ENj/Ysz-ov  are 
the  total  energies  of  the  Ni/YSZ  with  an  interface  oxygen  vacancy 
(Ni/YSZ-Ov)  with  and  without  the  CH  or  C  adsorbate,  respectively. 
The  Is2  of  H,  2s22p2  of  C,  2s22p4  of  O,  3d84s2  of  Ni,  4s24p65s24d2  of 
Zr  and  4s24p65s24d1  of  Y  are  treated  as  valence  electrons  in  the  DFT 
calculations.  The  Bader  charge  [54]  analysis  scheme  is  applied  to 
determine  the  atomic  charges  and  charge  transfer.  Pure  GGA 
functionals  underestimate  the  band  gap  and,  as  a  consequence, 
may  affect  the  other  properties.  This  drawback  can  be  remedied  in 
higher  levels  of  theory,  e.g.,  the  weighted  density  approximation 
(WDA),  screened  exchange  (sX)  [55],  GW  approximation  [56],  and 
hybrid  functionals  (HSE06,  PBE0,  B3LYP)  [57,58],  etc.  However, 
these  methods  demand  greater  computational  effort  and  are  not 
always  feasible  for  large  models  and  extensive  sampling.  For  this 
reason,  in  our  work  we  apply  the  GGA  approximation,  which  is 
known  to  give  good  energetics,  and  the  qualitative  description  of 
the  electronic  structure.  We  tried  to  use  the  international  system  of 
units  (SI)  in  this  paper.  However,  we  keep  some  of  the  popular  units 
used  the  micro  world  and  give  their  equivalent  in  SI  for  clarity,  e.g., 
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Fig.  1.  The  Ni/YSZ  model:  (a),  (b)  right  views,  and  (c)  side  views,  black  dash  square  in  (b)  represents  the  oxygen  vacancy  at  the  Ni/YSZ  interface,  black  dash  circles  in  (c)  represent  the 
intrinsic  oxygen  vacancies  in  the  YSZ  lattice. 


1  eV  =  1.60217733  x  10”19  J;  1  A  =  1CT10  m;  1 

e  =  1.60217733  x  10“19  C. 


3.  Results  and  discussions 

3.1.  The  adsorption  and  diffusion  of  CH  at  the  Ni/YSZ-Ov  interface 

Our  TPB  model  of  the  Ni/YSZ  system  contains  three  Ni  layers  on 
YSZ  as  shown  in  Fig.  1.  Since  oxygen  vacancy  usually  exists  after  the 
hydrogen  oxidation  at  the  Ni/YSZ  interface,  we  consider  the 
adsorption  and  diffusion  of  CH  at  the  Ni/YSZ-Ov  interface.  For  the 
Ni/YSZ-Ov  system,  there  is  a  significant  charge  transfer  of  0.96 
electrons  from  the  YSZ-Ov  support  to  the  Ni  part  as  compared  with 
that  in  Ni/YSZ  based  on  the  Bader  charge  analysis  [54  .  The  possible 
CH  adsorption  sites  are  shown  in  Fig.  1(c)  and  the  data  of  adsorption 
properties  are  shown  in  Table  1.  It  is  found  that  the  adsorption 
energy  and  excess  charges  of  CH  at  the  Ni  part  of  the  Ni/YSZ-Ov 
system  are  in  the  range  of  6.23-6.58  eV  and  0.53-0.66  e,  respec¬ 
tively,  which  are  close  to  the  values  (6.41  eV,  0.36  e)  of  CH  on  the 

Table  1 

The  adsorption  properties  of  CH  on  the  Ni/YSZ-Ov,  including  the  adsorption  energy 
(Fads,  in  eV),  the  bond  length  of  the  C-Ni(Zr)  (dC-Ni(Zr>  A)  and  C-H  (dC-H,  in  A),  the 
excess  number  of  electrons  of  the  adsorbed  CH,  CH-Chg  (in  e). 


Ni/YSZ-Ov 

Fads 

dc-Ni  (dc-Zr) 

dc-H 

CH-Chg 

1 

6.58 

1.82,  1.82,  1.90 

1.10 

0.67 

2 

6.41 

1.83,  1.83,  1.86 

1.10 

0.62 

3 

6.28 

1.81,  1.83,  1.91 

1.10 

0.62 

4 

6.38 

1.82,  1.82,  1.84 

1.10 

0.53 

5 

6.26 

1.82,  1.83,  1.84 

1.10 

0.58 

6 

6.23 

1.82,  1.84,  1.87 

1.10 

0.60 

7 

6.27 

1.82,  1.84,  1.85 

1.10 

0.61 

8 

6.30 

1.83,  1.83,  1.84 

1.10 

0.62 

Vac 

7.36 

1.91,  1.92  (2.36,  2.36) 

1.13 

1.23 

Ni(lll)  surface  [32].  In  contrast,  the  most  stable  CH  adsorption  is  at 
the  vacancy  site,  in  which  the  adsorption  energy  is  as  large  as 
7.36  eV  and  the  CH  gets  more  electrons  ( ~  1.23  e),  which  is  close  to 
the  value  of  the  sulfur  atom  at  the  vacancy  site  [42]. 

We  suppose  that  the  CH  diffusion  to  the  vacancy  site  experi¬ 
ences  a  two-step  process:  1)  the  CH  diffuses  on  the  ideal  Ni  (111) 
surface  beyond  the  TPB  region;  2)  the  CH  diffuses  from  the  Ni  part 
at  the  TPB  region  to  the  vacancy  site.  It  is  known  that  the  CH 
diffusion  on  the  ideal  Ni  (111 )  surface  is  quite  easy  (with  a  diffusion 
barrier  of  0.35  eV  [32]  from  the  HCP  site  to  the  FCC  site),  we  focus 
on  the  second  process  and  calculate  the  diffusion  barrier  of  CH  from 
the  most  stable  adsorption  site  1  to  the  vacancy  site. 


FS 


Fig.  2.  The  minimum  energy  path  (MEP)  for  the  diffusion  of  CH  fragment. 
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the  CH  fragment  can  diffuse  easily  to  the  vacancy  site.  The  large 
reverse  barrier  (0.91  eV)  means  that  the  CH  fragment  would  be 
trapped  at  the  vacancy  site. 

3.2.  The  dehydrogenation  of  CH  at  the  vacancy  site  of  the  Ni/YSZ-Ov 

Now,  we  study  the  dehydrogenation  of  CH  at  the  vacancy  site  of 
the  Ni/YSZ-Ov,  which  is  taken  as  the  initial  state  (IS).  The  converged 
dissociation  configurations  are  shown  in  Fig.  3.  It  is  found  the  H 
atom  is  not  favorable  at  the  O  top  site  at  the  interface  as  shown  in 
Fig.  3(a),  which  has  a  quite  large  endothermic  energy  (0.65  eV) 
compared  with  the  energy  of  CH  at  the  vacancy  site  of  the  TPB, 
while  the  H  atom  at  the  Ni  surface  is  relatively  more  favorable  when 
it  is  away  from  the  interface.  The  energies  of  configurations  shown 
in  Fig.  3(b)— (d)  are  (endothermic)  higher  in  energy  by  0.51,  0.12, 
0.09  eV  compared  with  the  energy  of  CH  at  the  vacancy  site.  The  CH 
dissociation  energy  profile  is  shown  in  Fig.  4  if  the  structures  in 
Fig.  3(a)— (d)  are  taken  as  the  final  states  (FS1,  FS2,  FS3,  FS4), 
respectively.  It  is  found  that  the  dissociation  of  CH  with  the  for¬ 
mation  of  OH  at  the  O  site  is  difficult  with  a  dissociation  barrier  of  as 
high  as  2.06  eV  (as  shown  in  the  transition  state  of  TS1),  while  the 
dissociation  paths  with  the  final  states  as  shown  in  Fig.  3(b)— (d)  all 
have  smaller  barriers  (0.74, 0.74, 0.83  eV)  as  shown  in  the  transition 
states  of  TS2,  TS3  and  TS4.  The  TS2  and  TS3  have  the  same  geometry 
configuration.  It  is  found  that  for  the  process  from  IS  to  FS3,  the 
dehydrogenation  of  CH  is  not  only  less  endothermic  with  smaller 
reaction  energy  (0.12  eV),  but  also  has  a  much  smaller  dissociation 
barrier  (0.74  eV)  as  compared  with  those  (0.54,  1.39  eV)  on  the 
Ni(lll)  surface  (e.g.,  the  region  beyond  the  TPB  region). 


Fig.  3.  The  final  states  for  the  CH  fragment  dissociation. 

The  minimum  energy  path  (MEP)  for  the  diffusion  of  CH  from 
site  1  to  the  vacancy  site  is  shown  in  Fig.  2.  The  diffusion  barrier  is 
only  0.13  eV  with  the  reverse  barrier  as  large  as  0.91  eV.  So  the  two 
processes  all  have  very  small  barriers  (0.35, 0.13  eV),  indicating  that 


4.  Discussion 

It  is  found  that  CH  would  easily  diffuse  to  the  vacancy  site  at  the 
Ni/YSZ-Ov  interface  and  the  dehydrogenation  barrier  (0.74  eV)  of 
the  CH  at  the  vacancy  site  is  lower  by  0.65  eV  as  compared  with  that 


Fig.  4.  The  MEPs  for  the  different  dissociation  paths  of  the  CH  fragment. 
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Fig.  5.  The  density  of  states  (DOS)  of  Ni  3d  at  the  Ni/YSZ-Ov  vacancy  site  and  on  the  Ni 
(111)  surface  (the  vertical  dash  line  represents  the  Fermi  level). 

(1.39  eV)  on  the  Ni  (111 )  surface.  The  density  of  states  (DOS)  of  Ni  3d 
at  vacancy  site  of  Ni/YSZ  and  that  on  the  Ni(lll )  surface  are  shown 
in  Fig.  5,  which  shows  that  the  Ni  3d  DOS  near  the  Fermi  level  is 
much  higher  at  the  Ni/YSZ  vacancy  site  than  that  on  the  Ni  (111) 
surface.  The  calculated  d-band  center  of  the  Ni  atom  at  the  vacancy 
site  is  much  shallower  than  that  on  the  Ni  (111)  surface  (-1.43 
vs  -1.66  eV).  The  lower  coordination  number  (6)  and  the  more 
excess  electrons  (0.46  e)  of  the  Ni  atom  at  the  vacancy  site  as 
compared  with  those  (9, 0.02  e)  on  the  Ni  (111 )  surface  make  the  Ni 
atom  at  the  vacancy  site  have  shallower  d  band-center  and  higher 
activity  for  the  dissociation  of  CFI.  The  high  active  Ni  atom  at  the 
vacancy  site  contributes  the  lower  dissociation  barrier  of  CH  (0.74, 
0.83  eV  vs  1.39  eV),  which  might  be  the  mechanism  for  the  carbon 
deposition  at  the  TPB  of  the  Ni/YSZ  anode. 

5.  Conclusion 

In  this  work,  the  carbon  deposition  at  the  TPB  of  the  Ni/YSZ-Ov 
system  is  studied  using  the  first-principles  method  based  on  den¬ 
sity  functional  theory.  It  is  found  that  the  CFI  fragment  (the  most 
stable  dissociation  product  of  CFI4  on  Ni  catalyst)  can  easily  diffuse 
and  be  trapped  at  interface  O  vacancy  site.  The  trapped  CH  can 
dissociate  to  C  and  H  atoms  with  lower  dissociation  barrier 
(0.74  eV)  as  compared  as  that  (1.39  eV)  on  the  pure  Ni  (111 )  surface. 
Therefore  the  interface  oxygen  vacancy  may  induce  the  formation 
of  carbon  deposition  at  the  TPB  of  the  Ni/YSZ  anode.  This  work 
sheds  new  lights  on  the  mechanism  on  the  carbon  deposition  of  the 
Ni/YSZ  anode. 
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